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ABSTRACT 

An  experimental  investigation  has  been  conducted  to 
examine  the  convective  heat  transfer  in  straight  and  curved 
ducts  of  rectangular  cross  section.   The  experimental  con- 
figuration was  modeled  as  infinite  parallel  plates  with  one 
wall  at  a  constant  heat  flux  and  the  opposite  wall  adia- 
batic.   Experiments  were  conducted   ;  steady  state  for 
turbulent  flow  using  average  Nusselt  numbers  to  compare  the 
heat  transfer  characteristics  of  the  straight  and  curved 
sections . 

The  development  of  Taylor-Gortler  vortices  in  the  curved 
section  proved  to  enhance  the  heat  transfer  rate  over  that 
of  the  straight  section.   Improved  heat  exchanger  and  turbine 
blade  design  could  both  result  from  a  better  understanding 
of  the  effects  of  curvature  on  the  rate  of  heat  transfer. 
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I.   INTRODUCTION 

A.   TAYLOR-GORTLER  VORTICES 

Since  the  early  part  of  the  twentieth  century,  consider- 
able research  has  shown  that  fully  developed  laminar  flow 
along  a  concave  wall  does  not  remain  two-dimensional  [Refs. 
1,  2,  3] .   The  flow  instead  forms  a  system  of  spiral  vortices, 
of  counter  rotating  pairs,  whose  axes  are  aligned  in  the 
principle  flow  direction.   This  phenomenon  is  the  result  of 
the  variations  in  the  centrifugal  forces  acting  on  the 
fluid  particles,  and  is  known  as  Taylor-Gortler  vortices. 
Figure  1  illustrates  the  type  of  fluid  motion  just 
described . 

In  a  channel  that  is  curved  in  the  streamwise  direction, 
those  fluid  particles  located  in  the  center  of  the  flow  cross- 
section  are  subjected  to  higher  centrifugal  forces  than  those 
fluid  particles  traveling  along  the  channel's  boundary  wall. 
As  a  result,  the  fluid  in  the  center  of  the  channel  moves  out- 
wardly towards  the  concave  boundary.   As  the  process  continues, 
the  fluid  particles  near  the  boundary  wall,  move  in  a  spanwise 
direction,  and  finally  radially  inward  replacing  the  outwardly 
moving  particles.   These  particles  then  come  under  the  same 
centrifugal  forces  and  the  process  repeats  itself.   The  re- 
sulting cyclic  motion  causes  the  formation  of  the  counter 
rotating  Taylor-Gortler  vortices,  considered  to  b  -  primai   / 

a  laminar  flow  phenomenon. 
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REPRODUCED    AT    GOVERNMENT    EXPENSE 


Primary    Flow,  Direction 


Figure  1.   Schematic  of  Taylor-Gortler  Vortices  in  a  Curved 
Channel . 
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It  has  been  observed  [Ref.  4]  that  the  heat  transfer  rate 
from  flow  along  a  concave  curved  wall  is  greater  than  that 
for  flow  along  a  straight  wall,  a  phenomenon  attributed  to 
the  additional  mixing  provided  by  the  secondary  motion  of 
the  Taylor-Gortler  vortices.   It  is  thought  that  the  cross- 
hatching  observed  on  reentry  vehicles  can  be  explained  in 
part  by  the  presence  of  streamwise  vortices  similar  to  the 
Taylor-Gortler  vortice  phenomenon. 

There  are  many  possible  applications  that  could  result 
from  a  more  thorough  understanding  of  Taylor-Gortler  vortices 
and  their  effect  on  heat  transfer  and  fluid  flow  character- 
istics.  Two  such  applications  could  lead  to  improved  heat 
exchanger  designs  and  improved  turbine  blade  cooling  [Refs. 
5,  6,  7] . 

B.   HISTORY 

The  instability  of  an  ir.   scid  fluid  flowing  past  a 
curved  boundary,  was  first  c..isidered  by  Lord  Rayleigh  _n 
1916  [Ref.  8].   By  assuming  that  the  fluid  was  non-viscous, 
he  determined  that  for  the  motion  of  an  inviscid  fluid  to 
remain  stable,  its  circulation  must  increase  with  increasing 
radius.   G.  I.  Taylor,  in  1923,  [Refs.  1,  9],  continued  this 
study  with  an  extensive  analytical  and  experimental  study 
of  viscous  fluids.   His  investigations  focused  on  the  flow 
between  two  cylinders,  in  which  the  inner  cylinder  rotated 
while  the  outer  cylinder  remained  stationary.   Taylor 
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ascertained  that  such  flows  become  unstable  when  the  value 
of  a  dimensionless  parameter  exceeded  a  critical  value  of 
41.3.   The  parameter,  known  as  the  Taylor  number  is  defined 
as : 


Ta  =  Re  / 


Ri 


where  ' d'  is  the  width  of  the  gap,  assumed  small  when  compared 
to  'Ri' ,  the  radius  of  the  inner  cylinder,  and  'Re'  is  the 
Reynolds  number  based  on  the  peripheral  velocity  of  the  inner 
cylinder.   Taylor  determined  that  for  those  cases  in  which 
the  value  of  the  Taylor  number  exceeded  the  critical  value, 
a  secondary  motion  developed  and  the  Taylor  vortices  formed. 
Figure  2  illustrates  this  fluid  motion. 

Instability  of  a  similar  nature  is  also  observed  when  a 
viscous  fluid  flows  in  a  curved  channel  due  to  a  pressure 
gradient  acting  along  the  channel  wall.   This  problem  was 
first  considered  analytically  by  W.  R.  Dean  [Ref.  10]  in  1928, 
for  a  channel  formed  by  two  concentric  cylinders,  where  the 
radius  of  the  inner  cylinder  was  large  in  comparison  to  the 
spacing  between  the  inner  and  outer  cylinder  walls.   Dean 
concluded  that  there  would  be  an  initiation  of  flow  instabil- 
ity, and  the  formation  of  vortices,  when  a  dimensionless 
parameter,  the  Dean  number,  exceeded  a  value  of  36.   The 
Dean  number  is  defined  as: 


De  =  Re  / 


Ri 
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where  ' d'  represents  the  channel  half-width,  ' Ri '  is  the 
inner  cylinder  radius,  and  'Re'  is  the  Reynolds  number  based 
on  the  mean  velocity  of  the  undisturbed  flow  and  the  channel 
half  width  'd'.   The  analytical  work  of  Dean  was  later  veri- 
fied by  W.  H.  Reid  [Ref.  11],  using  an  approximate  solution. 

In  1940,  H.  Gortler  [Ref.  2]  studied  the  stability  of 
laminar  boundary  layer  profiles  on  curved  walls  under  the 
influence  of  small  disturbances.   He  found  that  these  dis- 
turbances were  similar  to  the  vortices  studied  by  G.  I. 
Taylor.   Using  approximate  numerical  calculations,  Gortler 
concluded  that  the  disturbances  or  vortices  were  produced 
only  on  the  concave  boundary  walls  and  that  the  overall  flow 
profile  appeared  to  remain  laminar  in  nature.   These  results 
were  verified,  with  an  exact  solution,  by  G.  Hammerlin,  and 
reported  by  H.  Schlichting  [Ref.  12]  in  1955.   A.  M.  O.  Smith 
completed  an  even  more  extensive  numerical  analysis  that  fur- 
ther substantiated  these  findings  [Ref.  3].   The  results  of 
these  numerical  solutions  have  recently  been  demonstrated 
with  the  use  of  hot  wire  anemometry,  laser  doppler  systems, 
and  flow  visualization  techniques,  [Refs.  13,  14];  and  in 
1976,  Y.  Aihara  [Ref.  15]  conducted  a  non-linear  analysis 
of  Gortler  vortices. 

As  interest  began  to  develop  concerning  the  effects  of 
the  secondary  flows  associated  with  the  Taylor-Gortler 
vortices,  many  studies  were  published  concerning  the  in- 
fluences of  these  vortices  on  heat  transfer  in  the  laminar 
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and  turbulent  flow  regimes.   In  1955,  F.  Kreith  [Ref.  16], 
studied  th^  influence  of  heat  transfer  with  respect  to  the 
curvature  of  the  boundary  wall  for  fully  turbulent  flows. 
He  concluded  that  the  heat  transfer  from  the  heated  concave 
boundary  wall  was  considerably  higher  than  that  transferred 
from  the  convex  boundary  wall  of  the  same  curvature  and 
under  similar  turbulent  conditions. 

In  1965,  L.  Persen  [Ref.  17] ,  considered  the  ^oecial  cases 
of  very   igh  and  very  low  Prandtl  number  fluids  and  related 
the  inc    se  in  heat  transfer  from  a  curved  wall  to  the  pre- 
sence of  the  Taylor-Gortler  vortices. 

There  has  been  only  a  limited  amount  of  published  litera- 
ture dealing  with  the  flow  and  heat  transfer  in  curved  channels 
of  rectangular  cross-section  and  large  aspect  ratios,  where 
the  aspect  ratio  is  equal  to  the  spanwise  distance  of  the 
channel  divided  by  the  channel  height.   Much  of  what  has  been 
published  involves  the  development  of  numerical  approximations 
and  solutions  for  heat  and  mass  transfer  in  curved  ducts  of 
various  geometries.   K.  Cheng  and  M.  Akiyama  [Ref.  18]  devel- 
oped a  numerical  solution  for  forced  convection  heat  transfer 
with  laminar  flows  in  curved  channels  of  rectangular  cross- 
section,  but  only  for  small  aspect  ratios. 

In  1976,  A.  A.  Shibani  and  M.  N.  Ozisik  [Ref.  19],  using 
matched  asymptotic  expansion  techniques  for  a  wide  range  of 
Prandtl  n    ers,  solved  the  heat  transfer  probler,  between  paral- 
lel plates  tfith  turbulent  flow,  for  the  case  of  uniform  wall 

temperature. 

20 


Other  experimental  and  analytical  studies  that  are  worth 
citing  with  regard  to  this  present  study  follow.   Y.  Mori 

[Ref .  20] ,  obtained  results  for  hydrodynamically  fully  devel- 
oped flows  with  constant  wall  heat  flux,  in  curved  channels 
of  square  cross-section.   W.  M.  Kays  and  E.  Y.  Leung  [Ref.  21' 
reported  solutions  for  turbulent  flow  heat  transfer  in  a  con- 
centric circular  tube  annulus  with  a  fully  developed  velocity 
profile  and  constant  heat  rate  per  unit  length,  for  a  fluid 
of  Prandtl  number  0.7.   Results  for  large  aspect  ratio 
channels  with  rectangular  cross-sections  for  laminar  flows 
were  reported  by  M.  Durao  [Ref.  22]  and  J.  Ballard  [Ref.  23] . 
R.  Holihan,  Jr.  [Ref.  24]  reported  results  for  laminar  and 
transition  flows  and  S.  F.  Daughety  [Ref.  25]  reported 
results  for  turbulent  flows.   P.  F.  Brinich  and  R.  W.  Graham 

[Ref.  26]  reported  results  for  turbulent  flows  in  a  rectan- 
gular curved  channel  with  an  aspect  ratio  of  6,  for  the 
inner  wall  heated,  the  outer  wall  heated,  and  both  walls 
heated. 
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II.   INTENT  C   THE  STUDY 

The  purpose  of  this  investigation  was  to  examine  the 
effect  of  streamwise  curvature  on  the  heat  transfer  rate  in 
a  curved  rectangular  duct  of  large  aspect  ratio,  again  defin- 
ing aspect  ratio  as  the  spanwise  distance  divided  by  the 
channel  height.   The  flow  velocities  examined  were  in  the 
transitio:  and  turbulent  regime.   The  results  of  this  study 
were  compared  to  the  heat  transfer  rates  for  the  same  range 
of  velocities,  in  a  straight  duct  of  identical  aspect  ratio. 

Although  Taylor-Gortler  vortices  are  considered  to  be  a 
laminar  flow  phenomenon,  enhanced  heat  transfer  in  curved 
channels  for  transition  and  turbulent  flows  has  been  observed 
Even  though  it  has  not  been  proven  that  the  vortices  continue 
to  propagate  at  the  higher  flow  velocities,  this  enhancement 
in  the  transfer  of  heat  has  been  attributed  to  the  secondary 
flow  velocity  components  of  the  Taylor-Gortler   -ortices.   It 
is  believed  that  these  secondary  components  transport  the 
heated  fluid  from  the  outer  concave  wall,  inward  toward  the 
opposite  wall  of  the  channel,  displacing  the  cooler  fluid 
particles  and  causing  them  to  move  toward  the  heated  concave 
wall.   It  was  expected  that  similar  results  would  be  observed 
in  this  study. 

This  investigation  was  conducted  using  a  single  channel 
th  a  rectangular  cross-sectic.  and  constant  aspect  ratio. 

The  channel  incorporated  both  a  straight  test  section  and  a 
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curved  test  section.   The  results  obtained  in  each  test  section 
at  approximately  the  same  flow  rates  were  compared  in  an  effort 
to  determine  the  effects  of  the  Taylor-Gortler  vortices  on  the 
transfer  of  heat.   Also,  the  results  of  both  the  straight  and 
curved  sections  were  compared  to  the  results  of  Ballard  [Ref. 
23],  Holihan  [Ref.  24],  and  Daughety  [Ref.  25]. 

The  straight  section  results  of  this  study  were  compared 
to  the  Dittus-Boelter  equation  for  flows  in  circular  tubes 
[Ref.  26],  the  results  of  Kays  and  Leung  for  turbulent  and 
transition  flow  in  annular  passages  [Ref.  21] ,  and  the 
analytical  solutions  of  Petukhov  and  Popov  for  transition 
flows  in  circular  tubes  [Ref.  27] . 

The  curved  section  results  of  this  study  were  compared 
with  the  results  of  Brinich  and  Graham  [Ref.  28]  for  turbulent 
flow  in  a  rectangular  curved  channel,  and  also  to  the  solutions 
of  Petukhov  and  Popov,  mentioned  above. 
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III.   EXPERIMENTAL  WORK 

A.   DESCRIPTION  OF  THE  APPARATUS 

To  meet  the  objectives  of  this  investigation,  a  channel 
of  rectangular  cross-section  was  used.   Details  of  its  con- 
struction are  described  in  references  22  and  23.   Plexiglas, 
0.635  centimeters  thick,  formed  the  inner  and  outer  walls  of 
the  channel.   These  walls  were  separated  by  0.635  centimeter 
spacers  that  also  served  as  the  sides  of  the  channel. 

The  channel,  shown  in  Figure  3,  had  a  straight  section, 
122.0  centimeters  in  length,  followed  by  a  curved  section  of 
180  degrees  of  arc,  followed  by  another  straight  section,  91.4 
centimeters  in  length.   The  inner  curved  wall  of  the  channel 
had  a  radius  of  curvature  of  30.5  centimeters.   The  rectan- 
gular channel  was  0.635  centimeters  high  and  25.4  centimeters 
wide  resulting  in  an  aspect  ratio  of  40.   The  hydraulic  dia- 
meter was  1.229  centimeters.   The  cross-sectional  area  of 
the  channel  was  16.13  square  centimeters.   The  wetted  peri- 
meter was  52.07  centimeters.   Figure  4  shows  a  cross-sectional 
view  representing  both  the  straight  and  curved  sections  of 
the  channel. 

An  entrance  bell  constructed  of  plexiglas  was  connected 
to  the  straight  section  of  the  channel.   It  was  designed  and 
manufactured  according  to  ASME  nozzle  standards,  with  an 
elliptical  curved  base  on  a  major  axis  equal  to  ten  inches 
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and  a  minor  axis  of  one  inch.   Cheese-cloth  was  attached  to 
the  entrance  nozzle  to  prevent  foreign  matter  from  entering 
the  flow  channel.   An  aluminum  exhaust  nozzle  was  attached 
to  the  exit  of  the  channel  and  directed  the  flow  from  the 
channel  into  two  inch  diameter  pvc  piping.   The  two  inch  pvc 
piping  contained  a  measuring  section  which  could  be  fitted 
with  standard  ASME  concentric  orifices.   In  this  study  two 
orifice  plates  were  used.   An  orifice  with  a  diameter  of 
0.5325  inches  was  used  for  hydraulic  Reynold's  numbers 
below  5000,  and  an  orifice  with  a  diameter  of  1.065  inches 
was  used  for  hydraulic  Reynold's  numbers  above  5000.   Pressure 
taps  on  either  side  of  the  orifice,  at  1  and  1/2  pipe  diameters, 
and  connected  to  manometers  were  used  to  calculate  the  mass 
flow  rate  of  the  working  fluid,  which  was  air  at  room 
temperature . 

The  flow  of  air  was  drawn  through  the  channel  by  an 
electrically  driven  Spencer  Turbo  Compressor,  rated  30  horse- 
power at  3500  rpm,  and  550  cubic  feet  per  minute  at  70  degrees 
fahrenheit  and  one  atmosphere. 

The  channel  contained  two  test  sections  from  which  experi- 
mental heat  transfer  data  could  be  obtained.   The  straight 
test  section  was  located  120  centimeters  downstream  of  the 
entrance  bell  to  ensure  hydrodynamically  fully  developed 
flow.   The  straight  test  section  was  29.2  centimeters  in 
length  with  a  heated  test  section  area  of  741.7  square  centi- 
meters.  The  curved  test  section  was  located  in  the  lower 
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half  of  the  curved  portion  of  the  channel.   It  was  28.3 
centimeters  in  length,  subtending  an  arc  of  53.1  degrees, 
and  having  a  heated  test  section  area  of  718.4  square 
centimeters . 

In  each  of  the  heated  test  sections,  Temsheet,  a  carbon 
impregnated  porous  paper  with  a  uniform  electrical  resist- 
ivity, was  glued  to  the  outer  wall.   Joulean  heating  was 
used  to  heat  the  flow  of  air  through  the  channel.   Since  the 
electrical  resistance  of  the  Temsheet  is  not  constant,  but 
varies  slowly  with  temperature,  a  precision  resistor  with  an 
electrical  resistance  of  2.0262  ohms  was  connected  in  series 
with  the  Temsheet  to  allow  the  calculation  of  the  instantan- 
eous power  being  supplied.   Detailed  schematics  of  each  test 
section  are  shown  in  Figures  5  and  6. 

The  variables  that  were  measured  and  used  in  this  invest- 
igation were: 

(1)  the  temperature  of  the  air  entering  the  channel  (T.  ) 

(2)  the  temperature  of  the  air  leaving  the  channel  (T   ,  ) 

(3)  the  temperature  of  the  heated  outer  wall  for  each 

test  section  (T   ) 
wo 

(4)  the  temperature  of  the  unheated  inner  wall  for  each 

test  section  (T  . ) 
wi 

(5)  the  temperature  between  the  three  layers  of  insulation 

at  each  test  section  (T.   ) 

ins 

(6)  the  voltage  across  the  precision  resistor  (V  ) 

(7)  the  voltage  across  the  Temsheet  or  heater  (V,  ) 
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-Thermocouples 


Figure  6.   Curved  Test  Section,  Detailed  Schematic 
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(8)  the  temperature  of  the  air  at  the  orifice  (T   .^) 

orif 

(9)  the  pressure  upstream  of  the  orifice  (P   ) 

r  r  up 

(10)  the  difference  in  the  upstream  and  the  downstream 
pressure  across  the  orifice  (AP) 

(11)  the  atmospheric  pressure  (P    ) 

r       r  atm 

A  total  of  sixty-one  Copper-Constantan ,  30  gauge  thermo- 
couples were  located  throughout  the  test  apparatus  to  record 
the  desired  temperatures.   Readings  from  four  thermocouples 
spanning  the  width  of  the  channel  were  averaged  to  give  the 
temperature  at  the  entrance  and  exit  to  the  channel.   Five 
thermocouples  were  connected  in  parallel  and  inserted  between 
each  of  the  layers  of  insulation  at  each  of  the  test  sections. 
In  addition,  individual  thermocouples  were  positioned  in  each 
of  the  test  sections,  eight  in  direct  contact  with  the  Temsheet 
and  eight  in  the  plexiglas  of  the  unheated  inner  wall.   The 
beads  of  the  thermocouples  in  contact  with  the  Temsheet  were 
electrically  insulated  with  ENMAR  Heat  Resisting  Glyceryl 
Phthalate.   Small  diameter  holes  were  drilled  in  the  plexi- 
glas to  allow  the  installation  of  the  thermocouples.   One 
thermocouple  was  inserted  in  the  two  inch  pvc  pipe  downstream 
of  the  orifice  in  accordance  with  the  ASME  Power  Test  Code 
to  record  the  temperature  of  the  air  flowing  through  the 
orifice . 

The  thermal  insulation  consisted  of  1/2  inch  layers  of 
Armstrong  Armaflex  22  Sheet  Insulation,  a  flexible  foamed 
plastic  material.   By  positioning  thermocouples  between 
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these  insulation     ers,  the   eat  lost  to  the  environment 
could  be  computed   t  each  te~  -  section.   The  heat  loss 
through  insulation  surrounding  the  rest  of  the  channel  due 
to  the  heated  air  flowing  through  the  channel  was  small  in 
comparison.   See  Figure  7  for  a  detailed  sketch  of  the 
thermocouple  placement  and  attachment.   The  insulation  was 
held  in  place  by  the  use  of  velcro  strips  attached  to  the 
frame.   The  entire  channel  and  all  connections  between  the 
channel,  the  pvc  piping,  and  the  orifice  were  sealed  with 
General  Electrical  Silicone  Rubber  Sealant  Caulk  to  ensure 
there  was  no  leakage  of  air  into  the  channel  or  the  piping 
which  could  affect  the  recorded  temperatures  or  mass  flow 
rate. 

The  data  acquisition  system  used  for  this  experiment 
was  a  Hewlett  Packard  3054A  Automatic  Data  Acquisition/ 
Control  System  consisting  of  a  3456A  Digital  Voltmeter  and 
a  3497A  Data  Acquisition/Control  Unit.   Also  used  in  con- 
junction with  this  data  acquisition  system,  were  a  Hewlett 
Packard  9326  Computer  terminal  and  2671G  Printer.   Pressure 
measurements  were  taken  from  Meriam  vertical  manometers,  one 
whose  fluid  was  water  with  a  0  to  60  inch  range,  another 
using  water  with  a  0  to  3  0  inch  range,  and  a  third  whose 
fluid  was  mercury,  calibrated  to  read  inches  of  water  with 
a  range  of  0  to  115  inches.   A  photograph  of  the  channel 
and  associated    3t  equipment  are  shown  in  Figures  8,  9, 
and  10. 
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Sheets  of  cholesteric  liquid  crystal  slurries  with  a 
temperature  range  of  40-43°C  were  used  to  provide  qualita- 
tive measurements  of  the  surface  temperature  distribution 
on  the  Temsheet.   The  liquid  crystals  exhibit  dramatic  color 
changes  over  a  small  temperature  differential  and  allow  a 
reasonable  comparison  to  be  made  between  the  heat  transfer 
processes  on  each  test  section  by  observation  of  the  color 
pattern  displayed  for  varying  temperature  gradients  and  flow 
velocities.   Specifically,  the  occurrence  of  multi-colored 
stripes  aligned  in  the  flow  direction  on  the  curved  wall 
would  indicate  regions  of  high  and  low  temperatures  and  con- 
firm the  presence  of  Taylor-Gortler  vortices. 

The  addition  of  this  liquid  crystal  sheet  necessitated 
a  calculation  of  its  overall  heat  transfer  coefficient.   To 
determine  the  true  surface  temperature  of  the  heated  wall 
the  temperature  drop  across  the  sheet  must  be  known.   The 
sheet  was  composed  of  4.5  mils  of  a  liquid  crystal  slurry 
with  approximately  the  same  thermal  conductivity  as  body 
fat,  .5  mils  of  acrylic  backing,  and  a  2  mil  mylar  cover; 
resulting  in  an  overall  heat  transfer  coefficient  of 
522  w/°C.   Insulation  opposite  the  heated  sections  was 
fitted  with  removable  panels  to  allow  viewing  and  photo- 
graphing of  the  liquid  crystals. 

B.   EXPERIMENTAL  PROCEDURES 

The  experimental  procedures  followed  were  the  same  for 
the  straight  and  curved  test  sections.   The  results  obtained 
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from  the  straight  test  section  served  as  a  baseline  for  the 
comparison  of  the  curved  section  results.   Experiments  were 
conducted  by  varying  the  mass  flow  rate  of  air  corresponding 
to  Reynolds  numbers  from  approximately  800  to  25000,  where 
the  Reynolds  numbers  are  based  on  hydraulic  diameter.   For 
each  mass  flow  rate,  data  was  taken  automatically  by  the  data 
acquisition  system  and  the  data  was  immediately  reduced  and 
printed.   Preliminary  runs  were  performed  to  determine  the 
time  required  for  the  test  rig  to  come  to  steady  state. 
For  this  experiment  three  hours  were  allowed  for  steady 
state  to  be  reached.   After  three  hours,  data  was  taken  at 
ten  minute  intervals  to  ensure  that  a  steady  state  condition 
had  been  achieved.   The  criteria  for  steady  state  was  based 
on  three  variables: 

(1)  the  difference  between  inlet  and  outlet  temperature 

(2)  the  heated  boundary  wall  temperature 

(3)  the  unaccounted  heat  loss  in  the  system 

When  the  first  two  of  these  variables  varied  by  less  than 
two  percent  over  a  ten  minute  interval  and  the  third  began 
to  fluctuate,  it  was  considered  that  a  steady  state  condition 
had  been  achieved. 

It  was  also  necessary  to  determine  the  approximate  power 
setting  required  to  bring  the  heated  wall  of  each  test  section 
to  approximately  50  degrees  Centigrade.   It  had  been  deter- 
mined by  Holihan  [Ref.  24]  that  this  temperature  was  suffi- 
cient to  ensure  a  20  degree  Centigrade  difference  in  the 

heated  wall  and  the  unheated  wall. 
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For  each  run,  the  atmospheric  pressure,  the  pressure 
difference  across  the  orifice,  and  the  pressure  upstream  of 
the  orifice  were  entered  into  the  computer  for  the  mass  flow 
rate  calculations.   All  other  data  was  acquired  by  the  data 
acquisition  system  directly  from  the  experimental  apparatus. 

To  compute  the  instantaneous  power  supplied  the  following 
relationship  was  used: 

V   V 

Q    =     H    R 

^P      R 
F       PR 

The  precision  resistor  voltage  and  the  heater  voltage  were 
read  by  the  data  acquisition  system  and  the  resistance  of  the 
precision  resistor  was  a  known  constant.   A  diagram  of  this 
circuit  is  shown  in  Figure  11. 

Upon  completion  of  selected  data  runs  photographs  were 
taken  of  the  liquid  crystals  on  the  heated  section  to  provide 
a  visual  record  of  the  airflow  patterns. 

For  Reynolds  numbers  ranging  from  approximately  2400  to 
5600,  it  v/as  necessary  to  make  two  sets  of  data  runs  for  the 
straight  test  section.   The  first,  with  conditions  the  same 
as  for  other  flow  rates,  which  exhibited  a  pronounced  devia- 
tion from  the  near  linear  results;  and  a  second,  where  the 
flow  just  prior  to  the  straight  test  section  was  tripped 
with  sandpaper,  which  reestablished  the  near  linear  results. 
This  reproducible  anamoly  indicated  the  presence  of  a  tran- 
sition region  between  laminar  and  turbulent  flow. 
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Figure    11.       Circuit   Used    to   Measure    Power    Supplied    to 
the    Temsheet. 
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IV.   PRESENTATION  OF  DATA 

A.   ANALYSIS 

The  constant  heat  flux  surface  for  the  heated  wall  of  each 
test  section  was  provided  by  the  uniform  electrical  resis- 
tivity of  the  Temsheet.   The  insulated  unheated  wall  was  con- 
sidered adiabatic  since  the  heat  losses  through  that  wall  were 
negligible.   The  channel's  large  aspect  ratio  of  40  allowed 
the  modeling  of  the  channel  as  parallel  plates.   The  straight 
portion  of  the  channel  upstream  of  the  straight  test  section 
was  of  sufficient  length  to  ensure  that  the  flow  was  hydro- 
dynamically  developed  for  the  flow  velocities  of  this  study. 
The  straight  portion  of  the  channel  downstream  of  the  curved 
test  section  was  also  of  sufficient  length  to  ensure  that  the 
flow  exiting  the  channel  was  thoroughly  mixed  and  that  the 
temperature,  T    ,  was  an  average  bulk  temperature. 

Based  on  these  assumptions,  the  experimental  configura- 
tion was  modeled  as  forced  convection  between  parallel  plates, 
for  hydrodynamically  developed  and  thermally  developing  flows, 
subjected  to  a  constant  heat  flux.   The  boundary  conditions 
are  one  wall  at  a  constant  heat  flux  and  the  other  wall 
insulated. 

To  analyze  this  problem,  several  quantities  were  defined 
as  follows; 

The  heat  convected  to  the  air  was  calculated  using  the 

equation: 
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Q  .   =  m  C     (T   .  -  T.  ) 

air      £  Lr   out     in 


where  'C     '  was  the  specific  heat  of  the  air  at  constant 
pair  v 

pressure  and  *m'  was  the  mass  flow  rate  of  the  air.   Mass 
flow  rate  was  calculated  using  standard  ASME  flow  measure- 
ment techniques  [Ref.  33]. 

The  average  heat  transfer  coefficient  between  the  heated 
wall  and  the  flow  of  air  in  the  channel  was  def  .ed  by  the 

equation:  _ 

r-      air 
h   = 


ApL  AT 


where  '  Q  .  'is  defined  above.  '  A_.T  '  was  the  area  of  the 
air  '    PL 

Temsheet  in  the  test  section,  and  'AT'  was  the  difference 

between  the  arithmetic  average  heated  wall  temperature  (T   ) 

3  r  wo 

and  the  average  bulk  air  temperature  (T,,,).   The  average 

bulk  temperature  was  defined  as  the  arithmetic  mean  of  the 

entrance  and  exit  temperatures  (T.  ,  T   . ) . 

r  in'   out 

The  average  Nusselt  number  was  then  calculated  using: 

h  D,  , 
ST"  nd 

hd      K  . 
air 

In  this  equation  '  D,-,'  is  the  hydraulic  diameter  and  'K  .  ' 
n         hd  x  air 

is  the  thermal  conductivity  of  air. 

The  Reynolds  number  was  calculated  for  each  test  run  as 
follows : 

Re 


hd      Ay. 
c   air 


42 


where  again  'm'  and  ' D,  -,'    are  the  mass  flow  rate  and  hydrau- 
lic diameter  of  the  channel,  'A  '  is  the  cross-sectional 

c 

area  of  the  channel,  and  ' u    'is  the  dynamic  viscosity  of 

air 

the  air. 

For  the  curved  section  runs  the  Dean  number  was  defined  as: 


De  ■  Rehd  /  5T4 


where  ' Re,  , '  is  the  Reynolds  number  based  on  the  hydraulic 
diameter,  'D,  , '  is  the  hydraulic  diameter,  and  ' Ri '  is  the 
radius  of  curvature  of  the  unheated  inner  convex  wall. 

A  sketch  of  the  control  volume  and  a  set  of  sample  calcu- 
lations for  one  test  run  of  the  curved  section  are  given  in 
Appendix  B.   A  sample  calculation  for  the  uncertainty  analysis 
is  given  in  Appendix  A. 

B.   RESULTS 

The  data  obtained  from  each  experimental  run  was  evaluated 
utilizing  the  expressions  described  in  Part  A  above.   The 
major  parameters  resulting  from  this  evaluation  are  shown  in 
Tables  I  through  V.   Tables  I  through  III  contain  the  straight 
section  results  and  Tables  IV  and  V  contain  the  results  of 
the  curved  section.   A  plot  of  the  average  Nusselt  number 
verses  Reynolds  number  is  given  in  Figure  12  for  the  compari- 
son of  the  results.   Figure  12  shows  the  flow  passing  through 
the  transition  regime  of  the  straight  section  in  both  the 
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Table  I 

Summary  of  Straight  Test  Section  Results 
(Large  Orifice) 


Re,  , 

hd 

air 
(W) 

h 
(W/mC) 

AT 
(C) 

Nuhd 

4973 

47. 

09 

30. 

,80 

20. 

,61 

14, 

,97 

5633 

67. 

70 

35. 

.65 

25. 

,60 

17, 

.26 

6326 

83. 

83 

39. 

,19 

28. 

,83 

18, 

.81 

7096 

75. 

40 

44, 

.41 

22. 

,88 

21, 

.54 

7989 

98. 

82 

47. 

,99 

27. 

,75 

23, 

.09 

8935 

82. 

48 

52, 

,05 

21. 

,  36 

25, 

.10 

9969 

101. 

85 

56. 

,20 

24. 

,43 

27, 

.29 

11520 

121. 

29 

63. 

,12 

25. 

,90 

30, 

.52 

12772 

120. 

,00 

67. 

,01 

24. 

,13 

32, 

.35 

14357 

122. 

,20 

73, 

,52 

22. 

,40 

35, 

.43 

16137 

122. 

,09 

77, 

.78 

21. 

,16 

37. 

.50 

17990 

124. 

29 

83. 

,17 

20. 

,14 

40, 

.55 

20629 

125. 

,30 

89. 

,53 

18, 

,86 

43, 

.82 

24507 

127. 

,56 

101. 

,13 

17. 

,00 

48, 

.63 

Reynolds  number  and  Nusslet  number  are  based  on  hydraulic 
diameter . 
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Table  II 

Summary  of  Straight  Test  Section  Results 
(Small  Orifice  -  Untripped) 


Rehd 

Q  • 
air 

(W) 

h 
(W/mC) 

AT 
(C) 

hd 

1524 

27.17 

14.74 

24.84 

7.20 

2162 

32.19 

18.92 

22.93 

9.20 

2686 

34.81 

19.03 

24.65 

9.23 

3321 

39.19 

21.67 

24.37 

10.57 

4063 

40.31 

23.27 

23.34 

11.00 

4463 

39.28 

23.83 

22.21 

11.60 

4851 

52.09 

26.41 

26.58 

12.77 

5223 

51.15 

30.97 

22.25 

15.18 

Reynolds  number  and  Nusselt  number  are  based  on  hydraulic 
diameter. 
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1976 
2354 
2793 
3225 
3626 
4639 
5552 
7014 


Table  III 

Summary  of  Straight  Test  Section  Results 
(Small  Orifice  -  Tripped) 


Re,  ,  Q  .  h  AT  Nil,, 

hd  air        ,_7/  _.        ,_,x  Tia 

(W)        (w/mC)       (C) 


41.74 

17.91 

31.40 

8.62 

34.35 

19.87 

23.30 

9.62 

40.11 

22.55 

23.97 

10.88 

35.37 

24.52 

19.44 

11.95 

50.87 

28.29 

24.23 

13.72 

50.47 

33.84 

20.10 

16.44 

51.99 

39.33 

17.81 

19.18 

72.53 

45.85 

21.32 

22.27 

Reynolds  number  and  Nusselt  number  are  based  on    draulic 
diameter . 
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Table  IV 

Summary  of  Curved  Test  Section  Results 
(Large  Orifice) 


Rehd 

De 

Q  • 
air 

(W) 

E 
(W/mC) 

AT 
(C) 

Nuhd 

4952 

2257 

79. 

.49 

41. 

,71 

26. 

,11 

20. 

.39 

5650 

2574 

72. 

,36 

44. 

,95 

22. 

,05 

21, 

.97 

6560 

2989 

82. 

,47 

50. 

,05 

22. 

.57 

24. 

.20 

7578 

3453 

81. 

,40 

54. 

,48 

20. 

,47 

26. 

.39 

8444 

3847 

97. 

,95 

59. 

,04 

22. 

,73 

28. 

.56 

10000 

4555 

103. 

,57 

64. 

,62 

21. 

.96 

31, 

.11 

10972 

4999 

105. 

,79 

68. 

,81 

21. 

,06 

33, 

.80 

11515 

5256 

108. 

,32 

72. 

,74 

20. 

,40 

34. 

.97 

12067 

5498 

110. 

,84 

72. 

,50 

20. 

,94 

35. 

.41 

13000 

5922 

114. 

,52 

79. 

,05 

19. 

,84 

38. 

.19 

13723 

6252 

116. 

,24 

82. 

,86 

19. 

,22 

40, 

.34 

14493 

6602 

118. 

,59 

86. 

,79 

18. 

,72 

42. 

.36 

17238 

7854 

121. 

,95 

99. 

,76 

16. 

,75 

48, 

.41 

20517 

9348 

124. 

,65 

113. 

.84 

15. 

,00 

55, 

.50 

24208 

11029 

126. 

,31 

126. 

.43 

13. 

,68 

61. 

.20 

Reynolds  number  and  Nusselt  number  are  based  on  hydraulic 
diameter. 
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Table  V 

Summary  of  Curved  Test  Section  Results 
(Small  Orifice) 


Re,  , 
hd 

De 

Q  • 
air 

(W) 

h 
(W/mC) 

AT 
(C) 

hd 

1550 

706 

39. 

59 

22. 

75 

23. 

84 

10 

.86 

1828 

833 

40. 

55 

24. 

86 

22. 

35 

12 

.04 

2150 

978 

54. 

24 

27. 

41 

27. 

11 

13 

.10 

2426 

1105 

50. 

84 

29. 

66 

23. 

48 

14 

.27 

2751 

1254 

48. 

48 

30. 

38 

21. 

86 

14 

.83 

3372 

1536 

49. 

48 

33. 

56 

20. 

20 

16 

.00 

3815 

1738 

51. 

06 

36. 

39 

19. 

22 

17 

.65 

4312 

1964 

71. 

12 

39 

11 

24. 

91 

19 

.05 

5215 

2378 

72. 

52 

44 

02 

22. 

57 

21 

.22 

Reynolds  number  and  Nusselt  number  are  based  on  hydraulic 
diameter . 
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tripped  and  untripped  condition.   All  further  comparisons 
will  be  conducted  using  the  data  from  the  tripped  flow. 

The  results  indicate  an  increase  in  the  rate  of  heat 
transfer  with  increasing  Reynolds  number  for  both  the 
straight  and  curved  test  sections.   In  addition,  the  heat 
transfer  rate  was  higher  in  the  curved  section  for  each 
Reynolds  number  investigated. 

A  straight  line  least  squares  correlation  of  the  data 
for  the  present  study  resulted  in: 

Nu  =  0.02  Re,  ,°"79  2300  <  Re,,  <  7100 

hd  hd 

for  transition  flow  in  the  straight  test  section; 

Nu  =  0.22  Re,,0,53  1500  <  Re,,  <  5300 

hd  hd 

for  transition  flow  in  the  curved  test  section; 

Nu  =  0.07  Re,,0*65  8900  <  Re,  ,  <  25000 

hd  hd 

for  turbulent  flow  in  the  straight  test  section;  and 

Nu  =  0.03  Re,  -,0,76  8400  <  Re,  ,  <  25000 

hd  hd 

for  turbulent  flow  in  the  curved  test  section. 

In  the  short  laminar  regime  of  this  study  the  rate  of 
heat  transfer  was  about  forty-five  percent  higher  in  the 
curved  section,  indicating  the  presence  of  Taylor-Gortler 
vortices.   As  flow  progressed  through  the  transition  regime 
a  loss  of  enhancement  of  heat  transfer  rate,  from  forty  to 
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about  twenty  percent,  was  observed  as  flow  over  the  straight 
section  went  from  laminar  to  fully  turbulent.   In  the  low 
end  of  the  turbulent  regime,  Reynolds  number  between  54  00 
and  12000,  the  rate  of  heat  transfer  was  about  thirteen 
percent  higher  in  the  curved  section.   As  the  Reynolds 
number  continued  to  increase  beyond  12000  so  did  the  enhance- 
ment of  heat  transfer,  to  approximately  twenty-six  percent 
at  24000.   This  compares  favorably  with  earlier  studies  by 
Kreith  [Ref .  16] ,  where  he  reported  an  increase  in  the  rate 
of  heat  transfer  along  a  concave  wall  of  from  twenty-five 

to  sixty  percent  for  Reynolds  numbers,  based  on  hydraulic 

4        6 
diameter,  between  10   and  10".   Holihan  reported  an  increase 

of  fifteen  percent  for  laminar  flow  and  thirty  percent  for 

transition  flow  [Ref.  24].   Ballard  reported  an  increase  of 

eleven  percent  in  the  heat  transfer  rate  for  laminar  flows 

[Ref.  23] .   Daughety  reported  an  increase  of  twenty  percent 

in  the  heat  transfer  for  turbulent  flows  [Ref.  25] . 

The  actual  presence  of  Taylor-Gortler  vortices  was 

verified  by  the  observation  of  the  liquid  crystal  isotherm 

distribution  in  the  curved  section.   At  low  values  of 

Reynolds  number  and  through  most  of  the  transition  regime, 

liquid  crystal  bands  oriented  in  the  streamwise  direction 

were  visible  (Fig.  13).   As  the  Reynolds  number  increased 

turbulent  forces  began  to  dominate  in  the  curved  section, 

and  less  evidence  could  be  detected  for  the  existence  of 

the  vortices  (Fig.  14).   Over  the  entire  range  of  Reynolds 
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numbers  studied  for  the  straight  section,  the  liquid  crystals 
formed  essentially  straight  bands  perpendicular  to  the  flow 
direction,  indicating  approximately  two-dimensional  flow 
(Fig.  15).   Disruption  of  flow  in  the  transition  region 
due  to  the  insertion  of  the  sandpaper  is  illustrated  in 
Figures  16  and  17. 
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Figure  13.   Liquid  Crystal  Representation  of  Airflow, 
(Curved  Section,  Rehd  =  1100) 


Figure  14.   Liquid  Crystal  Representation  of  Airflow, 
(Curved  Section,  Rehd  =  9400) 
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Figure  15.   Liquid  Crystal  Representation  of  Airflow, 
(Straight  Section,  Re,d  =  12700) 
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Figure  16.   Liquid  Crystal  Representation  of  Airflow. 
(Straight  Section,  Rehd  =  4000,  Untripped) 


Figure  17.   Liquid  Crystal  Representation  of  Airflow, 
(Straight  Section,  Re.,  =  4000,  Tripped) 
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V.   DISCUSSION  AND  CONCLUSIONS 

Holihan  [Ref.  24],  in  his  study,  determined  that  there 
was  a  negligible  difference  in  the  fluid  bulk  temperature 
and  the  unheated  wall  temperature.   He  additionally  demon- 
strated that  the  radiated  heat  transfer  from  the  unheated 
wall  surface  was  minimal.   These  results  were  verified  by 
Daughety  [Ref.  25]  for  the  present  study  and  provide  the 
basis  for  the  assumption  that  the  heat  transfer  to  the  air 
flowing   through  the  channel  was  solely  by  convection  from 
the  heated  Temsheet. 

As  was  mentioned  earlier,  the  high  aspect  ratio  of  the 
channel  provided  the  basis  for  the  assumption  that  the  ex- 
perimental apparatus,  as  configured,  could  be  modeled  as 
infinite  parallel  plates.   Holihan's  data  [Ref.  24]  also 
tended  to  verify  this  assumption.   His  experimental  data, 
in  the  laminar  flow  region,  approached  the  theoretical  limit 
for  average  Nusselt  number  for  parallel  plates  with  one 
wall  heated  at  a  constant  heat  flux  and  the  opposite  wall 
adiabatic  [Ref.  29] . 

Based  on  the  assumptions  mentioned  above,  comparisons 
were  made  with  other  analytical  and  experimental  results 
that  were  of  the  same  or  similar  problem.   Again  the  problem 
being,  flow  between  infinite  parallel  plates  with  one  wall 
at  a  constant  heat  flux  and  the  opposite  wall  adiabatic. 
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A  comparison  of  the  experimental  results  of  this  study 
and  the  experimental  results  of  Ballard  [Ref.  23],  Holihan 
[Ref.  24]  and  Daughety  [Ref.  25]  are  shown  in  Figure  18, 
for  the  straight  section,  and  Figure  19  for  the  curved 
section.   Figure  20  shows  a  compilation  of  all  the  data  on 
one  plot.   Ballard's  data  was  strictly  in  the  laminar 
region;  Holihan *s  data  covered  the  laminar  as  well  as  tran- 
sition regions,  and  Daughety 's  data  was  strictly  in  the 
turbulent  region.   The  channel  aspect  ratio  was  40  for  each 
of  these  studies  and  the  experimental  procedures  were  similar 
in  each  case.   Ballard  and  Holihan  compared  their  data  to 
the  analytical  studies  by  McCuen,  [Ref.  30] ,  for  heat 
transfer  between  infinite  parallel  plates  with  constant 
wall  temperatures  and  heat  flux.   Daughety  compared  his 
data  with  the  results  of  Shibani  and  Ozisik  [Ref.  19] ,  and 
Kays  and  Leung  [Ref.  21].   Ballard's  and  Holihan's  data 
plotted  above  analytical  solutions  while  Daughety 's  data 
plotted  below.   Differences  were  attributed  to  the  variations 
in  the  geometries  used  in  the  studies,  side  wall  effects, 
inability  to  account  for  all  of  the  heat  transfer  processes, 
and  the  limitation  inherent  in  any  experimental  work. 

A  comparison  of  analytical  and  experimental  results 
with  those  of  the  present  study  for  the  straight  test  sec- 
tion is  shown  graphically  in  Figure  21.   The  results  of 
Petuhkov  and  Popov  [Ref.  27]  are  shown  for  turbulent  heat 
transfer  in  round  tubes.   The  correlation  used  was 
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Nu   =   (f/8)  Re  Pr 


(2/3) 


k1(f)+k2(Pr)/f/8(Prv"/  J,-l) 


where  'f  the  friction  factor,  has  been  calculated  using  the 
Filonenko  equation  given  in  Karlekar  and  Desmond  [Ref.  32]  as 

f  =  (1.82  log  Re  =  1.64)"2 


and 


kx(f)   =  1  +  1.4  f 
k„(Pr)  =  11.7  + 


Prd/3) 


The  Dittus-Boelter  equation  [Ref.  26]  for  heat  transfer  in 
a  straight  tube  for  constant  wall  temperature  is  plotted 
from: 

Nu  =  0.02  Re0*08 

For  a  Prandtl  number  equal  to  0.71.   Additionally,  the 
results  of  Kays  and  Leung  [Ref.  21]  for  heat  transfer  in 
annular  passages  are  shown  for  r*  =  1.0,  which  equates  to 
parallel  plates.   When  extrapolated  down  these  results 
continue  to  correlate  well  with  present  data. 

The  curved  section  data  is  plotted  in  Figure  22,  and 
is  compared  to  the  experimental  work  of  Brinich  and  Graham 
[Ref.  28] .   Also  plotted  as  a  reference  are  the  analytical 
solutions  of  Petukhov  and  Popov  [Ref.  27].   The  accuracy  of 
the  data  points  used  from  Brinich  and  Graham  are  subject  to 
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errors,  in  that  the  actual  values  are  not  given  in  their 
study  and  had  to  be  taken  from  a  plot  of  Stanton  number 
versus  Reynolds  number. 

Disregarding  the  small  differences  between  the  results 
and  correlations,  the  trend  appears  consistent.   Heat  trans- 
fer rates  can  be  seen  to  increase  both  with  increasing 
Reynolds  number  and  curvature  of  the  channel.   The  presence 
of  Taylor-Gortler  vortices  in  the  laminar  and  transition 
regime  seems  to  be  the  major  contributing  factor;  whereas 
increased  turbulence,  possibly  caused  by  the  formation  of 
these  vortices,  seems  to  dominate  at  higher  Reynolds  numbers. 
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VI.   RECOMMENDATIONS 

The  amount  of  experimental  work  remaining  in  the  study 
of  heat  transfer  in  curved  channels  is  significant.   The 
present  experimental  apparatus  has  been  used  over  a  range 
of  Reynolds  numbers  throughout  the  laminar,  transition  and 
into  the  turbulent  regime.   Additional  experiments  probing 
further  into  the  turbulent  regime  should  be  conducted  to 
better  understand  the  heat  transfer  processes  occurring 
there.   This  would  require  a  more  powerful  compressor  to 
generate  the  higher  Reynolds  numbers  and  a  larger  power 
supply  to  sufficiently  drive  the  heated  surfaces.   Finally, 
increasing  the  number  of  thermocouples  at  least  two-fold 
and  the  introduction  of  pressure  transducers  which  could 
be  connected  directly  to  the  data  acquisition  system  would 
both  prove  to  increase  the  accuracy  and  reliability  of  the 
experimental  apparatus. 

To  gain  further  insight  into  the  effect  of  curvature 
on  the  rate  of  heat  transfer,  additional  boundary  conditions 
for  the  test  should  also  be  considered.   An  apparatus  possess- 
ing a  heated  wall  on  each  side  of  the  test  section  might  be 
one  such  consideration. 
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APPENDIX  A:   EXPERIMENTAL  UNCERTAINTY 

The  uncertainties  for  the  major  variable  in  the  experi- 
ments, were  calculated  in  accordance  with  the  method  des- 
cribed by  S.  Kline  and  F.  McClintoch  [Ref.  29].   The  estimates 
of  the  uncertainty  in  the  measured  quantities  were  made  con- 
servatively.  As  a  result,  there  is  considerable  confidence 
in  the  uncertainties  as  calculated. 

The  following  equations  were  used  to  calculate  the 
uncertainties : 

2        2        2 

(1)  dm  =  m  [((  —  )       +  (£j£)   +  (S£)  ) 

A*  -i  9   1/2 

dp  •   2     -,  .  „  2 

+  (1/4  (<-SiS)  +    O    )] 

air 

,  1/2 

■A.    2     d  C   .   2     d(T    -T.  ) 

(0\    nn     -  n     i  tdm\  pair,     ,    out   in  >  . 

(2)  dQair  "  Qair  ( (1T)   +  (  C   .   )   +  (  T    -T . }  > 

pair         out   in 

1/2 
dQ  .   2     dAOT  2     ,.,  2    '  * 

(3)  dh  =  h  (<q-^)   +  (-^)   +  (^)  ) 

air  PL 

1/2 
Ar-   2  dD.,  2     dK  .   2  x/z 

(4)  <m  =  nu  ((^)    +  i-=p&)    +  (^ri£)   ) 

h         hd         air 

..  2     dDu,  2     dy     2     dA   2  1/2 
/rx  jt,       n     /  ,  dm.     ,   hd.     ,   air,     ,   c.  . 
(5)  dRehd  =  Rehd  (  (^   +  (-^— )   +  (-^— )   +  (— )  ) 

hd        air        c 
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The  values 

obtained 

for 

the  uncertainties  for  the  values 

in  the  curved 

section 

test 

run  at  a  Reynolds  number  of 

13000  are: 

Quantity 

Uncertainty 

A 

.0348 

APL 

.0019 

A 
c 

.0438 

C 

P 

.0020 

D,  , 
hd 

.0236 

h 

.0058 

K 

.0008 

K  . 
air 

.0038 

Nu 

.0127 

m 

.0050 

Q  • 
air 

.0058 

Rehd 

.0289 

T 

iblk 

.0004 

T. 
in 

.0004 

T     4- 

out 

.0004 

T 
wo 

.0002 

T   ,  -T. 
out   in 

.0027 

Y 

.0024 

AP 

.0500 

AT 

.0005 

air 

.0054 

P^v 

.0001 

air 
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The  major  source  of  uncertainty  in  the  average  Nusselt  number 
is  the  uncertainty  in  the  pressure  readings  from  the  mano- 
meter, which  affects  the  mass  flow  rate  calculation.   In 
order  to  obtain  a  more  accurate  pressure  measurement,  it 
has  been  recommended  that  pressure  transducers  be  used 
instead  of  vertical  manometers  in  future  studies. 
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REPRODUCED  AT  GOVERNMENT  EXPENSE 


APPENDIX  B:   SAMPLE  CALCULATIONS 

Figure  23  below,  shows  the  major  heat  transfer  components 
for  each  of  the  test  sections.   The  sample  calculations  that 
follow,  demonstrate  the  methods  used  by  the  computer  to  calcu- 
late these  components  as  well  as  the  Reynold's  number,  average 
heat  transfer  coefficient,  and  average  Nusselt  number  for  each 
set  of  data.   The  sample  calculations  are  for  the  curved  sec- 
tion, but  those  for  the  straight  section  are  similar.   Reynolds 
number,  Nusselt  number,  and  Dean  number  are  based  on  hydraulic 
diameter . 


°lo 


1 


2 


L -J--- 


Figure  23.   Energy  Balance  in  Straight  Section 
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A.   SAMPLE 

calculat: 

[ON  DATA 

VH 

= 

55.28 

V 

VR 

= 

4.53 

V 

RPR 

= 

2.0262  Q 

P   . 

atm 

= 

30.11 

in.  Hg 

AP 

= 

12.00 

in.  H20 

Pl 

= 

8.50 

in.  H20 

T. 

in 

= 

23.01 

C 

out 

= 

26.78 

C 

T 
wo 

= 

44.74 

C 

T  . 

Wl 

= 

24.64 

C 

insi 

= 

39.67 

C 

ins2 

= 

33.38 

C 

orif 

= 

26.57 

C 

K 
air 

= 

26.22 

x  10~3  W/mK 

C 
pair 

= 

1006  J/Kg  K 

y  ■ 
air 

= 

18.59 

x  10    Kg/ms 

K. 
ms 

= 

3.6495  x  10~2  W/mK 

AX. 
ins 

= 

0.01270  m 

e 
wo 

= 

0.70 

£   . 
Wl 

= 

0.40 

a 

= 

5.669  x  10"8  W/m2K4 

3 

= 

.5325 

Y 

= 

1.402 

gc 

= 

1  Kg 

m/NS2 

R 

= 

286.8 

Nm/Kg  K 
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F        =  1.0 
wo-wi 

ApL      =   .0730  m 

A        =   .0016  m 
c 

A  .      =   .002027  m2 
pipe 

D        =   .006  35  m 
c 

D  .      =   .0508  m 
pipe 

D   .  -    =   .027051  m 
orif 

P   ,     =   .5207  m 
wet 


B.   TEMPERATURE  CALCULATIONS 

1.   Bulk  Temperature  (T,,, ) 

blk 


Tbik  .   T"  *  T°"t  .  23.01  ±   26.78  .  24_9Q  c 


2.   Mean  Temperature  Difference  (AT) 


AT  =  T   -T,  ,,  =  44.74  -  24.90  =  19.85  C 
wo   blk 


3.   Temperature  Difference  in  the  Insulation  (AT.   ) 


AT.    =  T.   ,  -  T.   -  =  39.67  -  33.38  =  6.29  C 
ms     msl     ins2 


C.   POWER  CALCULATIONS 


1.   Power  Supplied  (Q  ) 
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O   =  V»  Vr  =  (55.28)  <4-53>  =  123  59  W 
^p     RpR        (2.0262)       '^-^  vv 


2.   Heat  Lost  Through  the  Outer  Plate  (Q,  ) 

3  lo 

AT  . 
ins _  6.29     

ylo    (AX.  J/(K.  J  (A„)  (0.0127)/ (3. 649  5x10"^)  (.0730 

ins     ins    xtJ-i 

=  1.32  W 


3.   Heat  Radiated  (Q  ) 

a.   Radiation  Resistance  (RD) 

1  -e         ,  1  -e  . 

R   -  w°  +      1       .  wi 

R    A„T  £       A^T  F      .  A„T  £  . 

PL  wo     PL   wo-wi  PL  wi 


[■—-  +  71-   "  1]  =  m  L,^[A  +  A-U=40.12  m-2 


A„T  L£      £  .     J    (0.0730)  L0.7  0.4 
PL    wo     wi 


b.   Heat  Radiated  (Q  ) 

c(T   2  -  T  .4) 
Q     _wo wi T[  =  ]K 

RR 

=  (5.669xl0~8)  [ (44.74  +  273)4  -  (  24 . 64  +  273 ) 4] 

40.12 


=  3.31  W 


D.   MASS  FLOW  RATE  CALCULATIONS 
1.   Pressure  Conversions 

a.   P    =30.11  in.  Hg  x  3374.1  =  101594.15  N/m2 

atm  3 
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b.  AP  =  12.0  in.  H20  x  248.64  =  2983.68  N/m2 

c.  P   =  (8.5  in.  HO  x  248.64)  +  101594.15  =  103707.59 

2 

N/m 

2.   Density  of  Air  (p    ) 

1  rair 


P-j  -    = 


Patm  (101594.15)  .     ,ooc    _     .    3 


air         R   T       .  , 
orif 


(286.8)  (26.57  +  273)         1.1825    Kg/nT 


3.   Expansion  Factor  (Y) 


Y  =  l-[0. 333+1. 145(62+0. 7B5+12313) ] 


YPX 


=  l-[ 0.33 3  +  1. 14 5  (  (.53  25)2  +  0.7( . 5325 ) 5  +  12 (. 5325 ) 13 ) ] 

2983.68 
x 


(1.402)  (103707.59) 
=  .9857 

4 .   Area  of  Orifice  (A) 

2 
.  "'Vif'  .  ^-02705)2  .  .0005747  m2 


5.   Mass  Flow  Rate  (m) 


m  =  YKA  /2gcp  . AP   =  (. 9857 ) K (. 0005747 ) 
3   air 


x  /2(1)  (1.1825)  (2983.68) 
m  =  .04759  K 
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Iterating : 

Assume  a  Reynolds  number       =  40000 
1  pipe 

Obtain  a  value  for  K,  the  flow  coefficient,  from 
reference  30.    K  =  .6350 

Solve  for  m.    m  =  .0302  Kg/s 

e  ,    *        n       m  Dpipe    (.0302)  (.0508) 

Solve  for  new  Re  .   =- r  r   =  -. — „„  on„_  .  .  ,  Q — =7; — =  __g 

pipe  A  .   y  .   (.002027)  (18.59x10  D 


pipe  air 


Re  .    =  40713 
pipe 


Check  convergence  and  repeat  process  if  necessary 
(Convergence  if  difference  less  than  .001) 
m  =  .0302  Kg/s 

E.   REYNOLDS  NUMBER  CALCULATIONS 


1    rp      =  ™  °Pipe    =  (.0302)   (.0508)      = 

x*     pipe    A  .   y  .     (.002027) (18.59xl0-b 
r   r     pipe  air 

9    po      ™  °c     .    (.0302)   (.00635) 

z*   Ked     A   y  .      (.0016) (18.59xl0-6~ 
c   air 

->  hd 

3 .   Re 


hd    Ay. 
c   air 

»    n       4  X  Ac  _  (4)  (.0016)  _ 

a.   Dhd  -    p      -  73207 '01229  m 

wet 


-   ^  Dhd    =  (.0302)  (.01229)    =  12478 

Kehd    A   y  .       (.0016) (18.59xl0-6 
c  Hair 


F.   HEAT  CONVECTED  TO  AIR  CALCULATION 

1.   Q  .   =  mC  (T    -T.  )  =  (.0302)  (1006)  (26. 78-23. 01) =114.  54  W 
wair      p   out   in 
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G.   AVERAGE  HEAT  TRANSFER  COEFFICIENT  CALCULATION 


7-     air  114.54 _Q  no  ..,  2^ 

h  =  ATT" AT"  =  (.0730X19.84)  =  79'°8  W/m  C 

C  i-l 


H.   AVERAGE  NUSSELT  NUMBER  CALCULATION 


m   =  ^  .  '79.08)  (.01229)  .  37_0? 

K  .       26.22  x  10 
air 


I .   DEAN  NUMBER  CALCULATION 


^     r,    /  Dc     .__.,_   /  .00635    co-7c 

De  =  Re   /  =■*-   =  40713  /  -nc   =  5875 

p    Ri  .305 
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APPENDIX  C:   CORRELATIONS 

The  correlations  obtained  for  the  present  study  were 
obtained  using  the  method  of  least  squares  for  a  first  degree 
polynomial,  as  outlined  by  C.  F.   Gerald  [Ref .  31] .   In  this 
method  the  values  of  Reynolds  number  (x)  and  Nusselt  number 

(Y)  were  first  converted  to  their  natural  logarithmic  value. 

2 

Next  the  summations  of  x . .  x .  ,  Y.,  and  x.Y.  were  calculated. 

1 '   l  '   l '       11 

The  values  of  these  quantities  were  then  placed  in  a  matrix: 


t  n  s%]  i  y  i  =  i  Yi  i 

Ex. Ex.       1       Ex.   Y. 
li  ii 


and  the  resulting  simultaneous  equations 


a  N  +  a, Ex.  =  EY. 
o      1   1      1 

a  Ex  .  +  a.  Ex  .   =  Ex  .  Y  . 
o   1      1   1        11 


were  solved  for  a   and  a,  . 

o       1 

The  equation  of  the  line  then  became: 


Ln  Nu   =  a   +  a,Ln  Re 
o     1 


or 


Nu   =  e  o    Re  1 


The  standard  error  for  the  present  study  correlation  was 
calculated  by  th   following  equation: 

76 


2       Z(Yi    V 

a   = 


N  -  n  -  1 


where:   Y.  is  the  actual  value  of  the  average  Nusselt  number 
obtained  experimentally. 

y.  is  the  value  of  the  average  Nusselt  number  calculated 
from  the  correlation  equation. 

N  is  the  number  of  data  points  used  in  the  correlation, 
n  is  the  degree  of  the  polynomial  used. 
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